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Abstract 

The insulin signaling patliway is critical for the control of blood glucose levels. Brown adipose tissue (BAT) has also been 
implicated as important in glucose homeostasis. The effect of short-term cold exposure on this pathway in BAT has not 
been explored. We evaluated the effect of 4 hours of cold exposure on the insulin pathway in the BAT of rats. Whole 
genomic microarray chips were used to examine the transcripts of the pathway in BAT of rats exposed to 4°C and 22 C for 
4 hours. The 4 most significantly altered pathways following 4 hours of cold exposure were the insulin signaling pathway, 
protein kinase A, PI3K/AKT and ERK/MAPK signaling. The insulin signaling pathway was the most affected. In the 
documented 142 genes of the insulin pathway, 42 transcripts (29.6%) responded significantly to this cold exposure with the 
least false discovery rate (Benjamini-Hochberg Multiple Testing: — loglO (p-value) =7.18). Twenty-seven genes (64%) were 
up-regulated, including the insulin receptor (Insr), insulin substrates 1 and 2 (Irsl and Irs2). Fifteen transcripts (36%) were 
down-regulated. Multiple transcripts of the primary target and secondary effector targets for the insulin signaling were also 
up-regulated, including those for carbohydrate metabolism. Using western blotting, we demonstrated that the cold 
induced higher Irs2, Irsl, and Akt-p protein levels in the BAT than in the BAT of controls maintained at room temperature, 
and higher Akt-p protein level in the muscle. Conclusion: this study demonstrated that 4 hours of cold exposure stimulated 
the insulin signaling pathway in the BAT and muscle of overnight fasted rats. This raises the possibility that acute cold 
stimulation may have potential to improve glucose clearance and insulin sensitivity. 
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Introduction 

Brown adipose tissue (BAT) is a thermogenic organ and 
consumes lipids and carbohydrates. BAT metabolic activity is 
increased in a cold environment [1,2]. The unique energy 
consumption of stimulated BAT might be useful in controUing 
obesity and diabetes by modification of the body's energy balance 
and lead to more expenditure of energy and less deposition of fat. 
Over a decade ago functioning BAT was unexpectedly identified 
in adult humans by fluodeoxyglucose positron emission tomogra- 
phy ( F-FDG PET), a modern functional imaging modality 
[3,4,5] . Adults with FDG-avid BAT have lower body mass indices 
(BMI) than individuals with non-FDG avid BAT [3,6]. Adults 
without FDG-avid brown adipose tissue on PET imaging had a 
higher risk of abnormally increased glucose levels than patients 
with FDG-avid brown adipose tissue [7] . This finding leads to the 
hypothesis that purposely stimulating BAT could have a role in 
controlling obesity and diabetes. 

Cold temperature is a natural stimulator of BAT thermogenesis 
[4]. Cypess et al reported that cold activated human brown 
adipose tissue but sympathomimetics (ephedrine) did not [8] . In a 
series of rat experiments, cold exposure (48 hours at 4° or 5°C) 
reproducibly activated brown adipose tissue, and was associated 
with accelerated glucose clearance [9,10,11,12]. 



Gasparetti et al addressed the correlation between glucose 
clearance and the insulin signaling pathway in brown and white 
adipose tissue, skeletal muscle, and liver of rats chronically exposed 
to a cold environment (at 4°C for 8 days) or at room temperature 
(23°C) [13]. They focused on initial and intermediate steps of the 
insulin signaling pathway using immunoblotting and immunopre- 
cipitation. AU of the tissues except the liver increased glucose 
uptake after the cold exposure. Only brown adipose tissue 
responded to the cold exposure by significandy upregulating the 
protein level of insulin receptor substrate 2 (IRS2). This cold 
condition significantly potentiated insulin-induced phosphoryla- 
tion of the insulin receptor, insulin receptor substrate 1 and insulin 
receptor substrate 2 in the BAT of the rats. The cold also increased 
glucose clearance during a glucose tolerance test. These findings 
support increased insulin sensitivity in the BAT. The cold-induced 
modulation of insulin signaling caused BAT to use energy for 
thermogenesis and led to an overall increase in glucose utilization 
[13]. 

In previous rat experiments, we reproducibly demonstrated that 
acute cold exposure (4°C for 4 hours) significandy produced 
higher FDG uptake in the BAT than the other tissues of rats. The 
FDG uptake was also much higher than that found in the non- 
stimulated BAT of the rats at room temperature [14]. These data 
led to our hypothesis that only 4-hours of cold exposure might 
induce the changes in the transcripts involved in the insulin 
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receptor signaling pathway of BAT. We assumed that the more 
transcripts in the pathway that are significantly regulated, the 
more reliable the prediction of the status of the pathway would be. 
Since no .study has focused on all of the transcripts involved in the 
insulin signaling pathway in the BAT of rats, we used updated rat 
whole genomic chips to assess the responses of all transcripts 
involved in the insulin signaling pathway and the associated 
carbohydrate metabolism in BAT of rats after 4 hours of cold 
exposure at 4°C. 

Materials and Methods 

Animals 

Female Lewis rats (body weight 126-150 g) were purchased 
from Charles River Laboratories, Inc (Wilmington, MA). They 
were housed in the animal facilities (room temperature 22-23°C) 
and exposed to a 1 2 hour light/ dark cycle at least 1 week before 
the start of experimentation. During this adaptation period, food 
and water were given ad libitum. The study was performed as part 
of a project approved by the animal research committee at the 
Johns Hopkins Medical Institutions. 

BAT activation 

Female rats (n = 6) were randomly divided into two groups of 
three each. The rats in the control group (CTL) were kept at room 
temperature, 22°C. The rats in the treatment group were placed at 
4°C for 4 hours [14,15,16]. In both groups, food was withdrawn 
for the 4-hour conditioning period. Animals were humanely 
sacrificed using CO2 asph^'xiation. At sacrifice, the BAT in 4°C 
cold-treated group appeared highly perfused and was darker and 
redder in color than the BAT from the 22°C control group, 
consistent with BAT activation after 4°C exposure. 

RNA extraction and microarray-based assessment of 
relative transcript abundance 

Interscapular BAT was dissected from each animal immediately 
after the animals were sacrificed and separately stored in liquid 

nitrogen. Total RNA was extracted from the samples with the 
RNeasy kit purchased from QIAGEN (German town, MD). The 
quality and cjuantit)' of the RNA were evaluated with a NanoDrop 
UV-Vis spectrophotometer and Agilent Bioanalyzer. The Rat 
Expression 230 2.0 microarray chip from Afiymetrix (Santa Clara, 
CA) is a whole-genome array to interrogate 31,099 transcripts and 
variants from the rat genome, and utilized by the Microarray Core 
Laboratory' at Johns Hopkins Medical Institutes. A total of 6 chips 
were used (3/group), one per rat BAT sample. 

Analysis of the microarray data and signaling pathways 

To remove sources of variation of non-biological origin among 
arrays, the microarray raw data were first normalized by 
RMAexpress (Robust Multichip Average express), software [17]. 
Tli(" significance of gene expression levels between the cold and 
control groups were assessed by Spotfire (TIBCO, Somerville, 
MA), commercial software for gene array analysis. 

The transcripts with a p-value less than 0.01 were selected for 
further analysis of signaling pathways. Using analytic software. 
Ingenuity Pathway Analysis (IPA, Ingenuity Systems, Redwood 
City, CA), th(-se transcripts were grouped based on their biological 
functions and apparent involvement in signaling pathways. The 
Fisher's Exact Test was used to compute the p-value with the 
Benjamini-Hochberg multiple testing correction applied to 
control for the false discovery rate. Since most p-values from 
the Fisher's and Benjamini tests were very small, we present them 
as — logio (p-value), which we term the "p-score". 



To assess the activation states of various biological processes, 
activation z-scores were computed based on the change in the 
direction of the transcripts and the interactive relationship among 
genes. If the changes in our data were consistent with the changes 
of the known activation model in the Ingenuity System, the z- 
scores would be positive. The z-scores would be negative if 
opposite occurred. In the current study, z-scores would be 
considered significant if they were greater than 2 (activation) or 
smaller than -2 (inhibition). The reports of statistical analysis, 
pathways, and biological processes were generated by the use of 
IPA (Ingenuity Systems, www.ingenuity.com, Redwood City, 
CA). 

Western blotting for insulin signaling 

Nine rats were fasted overnight. The next morning, they 
were randomly put into one of three conditions: room 
temperature -I- saline (RT, nonstimulated), cold temperature -I- 
saline (4°C for 4 hours), and room temperature with insulin 
administration (0.75 units/kg), HumulinR from Eli Lilly, 
USA[18], as positive controls of insulin signaling activation. 
All injections (saline or insulin) were intraperitoneally admin- 
istrated 10 minutes before sacrifice. The intrascapular brown 
adipose tissue and scapular muscles were immediately dissected 
and separately frozen with dry ice. At the end of the 
experimental condition each rat was humanely sacrificed with 
the method described above. 

Approximately 10 mg brown adipose tissue was used for the 
preparation of lysate from each sample with RIPA lysis buffer 
(9806s, Cell Signaling, MA, USA) plus phosphatase inhibitor. 
Muscle lysates were prepared following the same steps. 

Approximately 30 |J,g protein lysate per lane was loaded to 
NuPAGE Bis-Tris Mini Gels (MAN0003679, life tech, CA, USA). 
After electrophresis, the samples were blotted onto nitrocellulose 
membranes. The membranes were visualized by Ponceaus 
staining (P3504, Sigma, MO, USA) and cut into distinct strips 
based on the target protein's molecular weight. Rabbit primary 
antibodies against the insulin receptor (#3025), Irs2 (#4502), Irsl 
(#2382), Akt-p (#9271), Akt-t (#9272) were obtained from CeU 
Signahng. Anti-Glut4 (SAB4300667), anti-P actin (A1978, 
Sigma), and actin (A2103, Sigma) were obtained from Sigma. 
Donkey anti-rabbit Ig IIRP conjugate (NA9340) was obtained 
from Amersham (GE Healthcare, Buckinghamshire, UK). Signals 
were detected by enhanced chemUuminescence plus the Western 
Blotting Detection System (RPN2232, GE Healthcare, Buck- 
inghamshire, U.K.). 

After an overnight incubation with block buffer (TBST+5% 
nonfat), the primary antibodies 1:1000 except for fi-actin (1:5000) 
and actin (1:35000) were incubated with the blotted nitrocellulose 
membranes overnight. These processes were performed at 4°C 
with gentle agitation on an orbital shaker. The secondary 
antibodies (1:2000 for anti-rabbit-HRP and 1:5000 for anti- 
mouse-HRP) were incubated for 1 hour at room temperature with 
gentle agitation. The signals on Hyperfilm ECL (#28906835, GE 
healthcare, Buckinghamshire, UK) were scanned into a computer 
and analyzed with Image Studio Lite Ver 3. 1 (LI-COR Corporate, 
NE, USA). The fi-actin signal in each sample served as a 
denominator to normalize the insulin signaling in the current 
experiments. We obtained relative signal intensities for further 
comparisons to determine fold changes from controls. Student's t 
test was employed to examine the protein expression differences 
among groups. 
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Results 

Cold-induced change in the insulin receptor signaling 
pathway 

The top 4 altered canonical pathways in BAT responding to 
cold stimulation were insulin receptor signaling, protein kinase A 
signaling, PI3K/AKT signaling, and ERK/MAPK signaling. AU 
of them had low false discovery rates (p<0.001). 

In the insulin receptor signaling pathway with the documented 
142 genes, 42 transcripts (29.6%) significantly responded to 
4 hours of cold exposure. Among 42 genes, the down-regxilation 
was 15/42 (36%) and the up-regulation was almost double, 27/42 
(64%). 

To visualize the relationship among the transcripts, these 42 
transcripts were graphed to a modified diagram of the insulin 
receptor signaling pathway with red (up-regulation), green (down- 
regulation), or non-color (unchanged) by IPA (Figure 1). The 
diagram revealed concurrent up-regulation in the insulin receptor- 
MAPK branch, but mixed regulation in the insulin receptor-PI3K 
branch (Figure 1). 



At the initial part of the pathway, the transcript of the insulin 
receptor (Insr) was up-regulated, and associated with concurrendy 
up-regulated transcripts of insulin receptor substrate 1 (Irsl), 
insulin receptor substrate 2 (Irs2), and Jakl (Janus kinase 1) 
(Table 1 and Figure 1). The transcripts of phosphatases, which 
dephosphorylate the insulin receptor, were not changed. 

In the MAPK branch, concurrent up-regulation is shown 
(Figwe 2), including Sosl (son of sevenless homolog 1), Nek (non- 
catalytic region of tyrosine kinase adaptor protein), Crk (V-Crk 
sarcoma virus CTIO), Ras (rat sarcoma), Raf (rapidly accelerated 
fibrosarcoma), growth factor receptor bound protein 10 (GrblO), 
Fyn (FYN oncogene related to Src, Fgr, Yes) (Table 1 and 
Figure 1). The transcript of MAPK8 (mitogen-activated protein 
kinase 8)/Jnkl was up-regulated. 

We analyzed the three classes of the PI3K family, class I, II, and 
III. The Pik3ca (catalytic, class I) and Pik3c2a (class II) were up- 
regulated. The Pik3r2 (regulatory subunit 2, class I) and Pik3c3 
(class III) were down-regulated. However, the sum of phospho- 
inositide-3-kinase (Pik3) was down-regulated since the down- 
regulation of Pik3r2 and Pik3c3 was greater than the up-regulation 
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Figure 1 . Diagram of the insulin signaling pathway modified by coding the change of transcripts in the BAT of rats responding to 4- 
hours of cold stimulation (4°C) as compared to the control group (22°C). Up-regulated transcripts are red (the darker red, the more up- 
regulation). Down-regulated ones are green (the darker green, the more down-regulation). The unchanged transcripts have no color. The signs 
between two components represent activation (arrow) and inhibition (bar or bar crossing arrow). Our data show that 4-hours of cold stimulation 
regulated multiple transcripts of the insulin signaling pathway in BAT of rats. 
doi:1 0.1 371/journal.pone.0099772.g001 
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Table 1. Insulin receptor signaling pathway. 





Symbol 


Gene Name 


Log Ratio 


p-value 


Type(s) 


Affymetrix 


SGKl 


serum/glucocorticoid regulated kinase 1 


3.34 


6.03 E-05 


kinase 


1367802_at 


IRS2 


insulin receptor substrate 2 


2.40 


3.27E-05 


enzyme 


1371091_at 


INPP5B 


inositol polyphosphate-5-phosphatase, 75kDa 


2.06 


6.88E-04 


phosphatase 


1 388502_at 


FYN 


FYN oncogene related to SRC, FGR, YES 


1.51 


6.65 E-05 


kinase 


1373683_at 


IRS1 


insulin receptor substrate 1 


1.42 


7.73 E-04 


enzyme 


1369771_at 


INSR 


insulin receptor 


1.38 


2.96E-03 


kinase 


1 392043_at 


KRAS 


v-Ki-ras2 Kirsten rat sarcoma viral oncogene homolog 


1.15 


7.65 E-06 


enzyme 


1370035_at 


CRK 


v-crk sarcoma virus CT10 oncogene homolog (avian) 


1.07 


9.14E-05 


other 


1 392462_at 


PRKARIA 


protein kinase, cAMP-dependent, regulatory, type 1, 
alpha 


1.04 


5.22E-04 


kinase 


1 386905_at 


PIK3C2A 


phosphatldylinosltol-4-phosphate 3-kinase, catalytic 
subunit type 2 alpha 


0.87 


4.36E-04 


kinase 


1 379433_at 


PTEN 


phosphatase and tensin homolog 


0.80 


1 .30E-03 


phosphatase 


1370112_at 


SH2B2 


SH2B adaptor protein 2 


0.78 


2.00E-03 


other 


1 368605_at 


RRAS2 


related RAS viral (r-ras) oncogene homolog 2 


0.74 


2.82E-03 


enzyme 


1 382058_at 


PDPKl 


3-phospholnositlde dependent protein kinase-1 


0.72 


5.38E-03 


kinase 


1370052_at 


NCKl 


NCK adaptor protein 1 


0.68 


2.87 E-03 


kinase 


1373940_at 


PPP1R3D 


protein phosphatase 1, regulatory subunit 3D 


0.66 


7.79E-03 


phosphatase 


1373656_at 


GRB2 


growth factor receptor-bound protein 2 


0.53 


5.95E-03 


other 


1368385_a_at 


GRB10 


growth factor receptor-bound protein 10 


0.51 


6.56E-03 


other 


1371517_at 


PPP1R12A 


protein phosphatase 1, regulatory subunit 12A 


0.51 


8.78E-03 


phosphatase 


1 382307_at 


MAPK8 


mitogen-activated protein kinase 8 


0.49 


1.27E-03 


kinase 


1 37961 2_at 


SOSl 


son of sevenless homolog 1 (Drosophila) 


0.48 


3.52E-03 


other 


1389710_at 


EIF4E 


eukaryotic translation initiation factor 4E 


0.47 


5.09E-03 


translation regulator 


1 398799_at 


SYNJl 


synaptojanin 1 


0.43 


1 .35E-03 


phosphatase 


1 370070_at 


PDE3B 


phosphodiesterase 3B, cGMP-inhibited 


0.37 


5.71 E-03 


enzyme 


1369157_at 


RAFl 


v-raf-1 murine leukemia viral oncogene homolog 1 


0.34 


1.37E-03 


kinase 


1369932_a_at 


PIK3CA 


phosphatidylinosltol-4,5-blsphosphate 3-kinase, 
catalytic subunit alpha 


0.27 


9.47E-03 


kinase 


1 382366_at 


JAKl 


Janus kinase 1 


0.26 


1.62 E-03 


kinase 


1 383478_at 


STX4 


syntaxin 4 


-0.16 


4.38E-03 


transporter 


1370014_at 


F0X04 


forkhead box 04 


-0.26 


1.47 E-03 


transcription regulator 


1 372652_at 


EIF2B4 


eukaryotic translation initiation factor 2B, subunit 4 
delta, 67kDa 


-0.28 


8.86E-03 


translation regulator 


1 386970_at 


PRKACA 


protein kinase, cAMP-dependent, catalytic, alpha 


-0.30 


7.88E-04 


kinase 


1371578_at 


PIK3C3 


phosphatldylinositol 3-klnase, catalytic subunit type 3 


-0.40 


1.17E-03 


kinase 


1369655_at 


INPPLl 


inositol polyphosphate phosphatase-like 1 


-0.41 


7.55E-03 


phosphatase 


1 368022_at 


EIF4EBP1 


eukaryotic translation initiation factor 4E binding 
protein 1 


-0.43 


6.92E-03 


translation regulator 


1 386888_at 


PPP1R3C 


protein phosphatase 1, regulatory subunit 3C 


-0.48 


1.93 E-03 


phosphatase 


1373108_at 


AKT2 


v-akt murine thymoma viral oncogene homolog 2 


-0.61 


1.31 E-03 


kinase 


1 388765_at 


PKKAG2 


protein kinase, AMP-actlvated, gamma 2 non-catalytic 
subunit 


-0.62 


1 .38E-03 


kinase 


1 373952_at 


INPP5K 


inositol polyphosphate-5-phosphatase K 


-0.64 


1.12E-04 


phosphatase 


1373627_at 


PRKAR2B 


protein kinase, cAMP-dependent, regulatory, type II, 
beta 


-0.68 


1 .53E-04 


kinase 


1371133_a_at 


PIK3R2 


phosphoinosltide-3-klnase, regulatory subunit 2 (beta) 


-0.71 


2.44E-04 


kinase 


1376190_at 


LIPE 


lipase, hormone-sensitive 


-0.73 


6.29E-03 


enzyme 


1387132_at 


INPP5F 


inositol polyphosphate-5-phosphatase F 


-0.85 


7.41 E-03 


phosphatase 


1389176_at 



The cold-Induced transcripts (p<0.01) of the Insulin signaling pathway are shown above. The interactive relationships among the transcripts are shown in Figure 1. 
Based on the extent of the changes, the transcripts are listed from the highest to the lowest. Log Ratio = log2 (average cold/average control). The p-values are 
represented in the scientific format of Microsoft Excel (for example, 6.03E-05 = 6.03x10"°^). 
doi:l 0.1 371 /journal.pone.0099772.t001 
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of Pik3ca and Pik3c3 (Table 1 and Figure 2). Pdkl (3- 
phosphoinositide dependent protein kinase- 1), a downstream 
kinase of PI3K, was up-regulated, associated Sgkl (serum/ 
glucocorticoid regulated kinase 1) up-regulation but the Pkc family 
was unchanged. Akt2 (v-akt murine thymoma viral oncogene 
homolog 2) was down-regulated in the current study. Transcripts 
of Pka and hormone sensitive lipase (Lipe) were down-regulated in 
the setting of up-regulated PdeSb (phosphodiesterase 3B, cGMP- 
inhibited). Overall, the over-expressed transcripts in the PI3K/ 
AKT branch were greater than the underexpressed transcripts. 

Cold-induced gene change involved in carbohydrate 
metabolism 

Using IPA Downstream Effects Analysis, we calculated p-values 
and z-scores to predict the association among the changed 
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Figure 2. Western blot showing expression of the insulin 
receptor, Irs2, Irs1, Akt-p, Akt-t, Glut4, and Actin in the BAT of 
rats with/without stimulation. Nine rats were fasted overnight. The 
next morning, they were randomly exposed to one of three conditions: 
room temperature + saline (RT, control), cold temperature + saline (4°C 
for 4 hours. Cold), and room temperature with insulin administration 
(0.75 units/kg, RT-l-ins). The preparation of BAT lysates was performed 
as described in Materials and Methods. Panel A: Protein levels of Insulin 
receptor, Irs2, Irsl, Akt-p, Akt-t, Glut4, and Actin, analyzed and 
compared to controls. Panel B: Columns represent mean ± SEM 
protein levels determined from "A" as a percentage control values of 
three experiments (t-test with a = 0.05, * p<0.05). The insulin signaling 
pathway in the BAT of rats post the acute cold exposure was more 
active than controls with increased Irs2, Irsl, Akt-p post cold exposure. 
doi:1 0.1 371/journal.pone.0099772.g002 



transcripts in our data as compared to the transcript interaction 
related to carbohydrate metabolism (Table 2). There were two sets 
of transcripts corresponding to carbohydrate metabolism. One set 
of transcripts is involved in the metabolism of carbohydrates (p- 
value 5.16x10 and had an activated z-score 2.643 giving 66% 
up-regulation and 34% down-regulation in 166 transcripts). The 
other set of transcripts is related to synthesis of carbohydrate (p- 
value 7.60x10 and an activated z-score of 2.556 due to 70% 
up-regulation and 30% down-regulation in 125 transcripts). 

Cold-induced change in glucose transporters 

Glucose transporters within plasma membranes are crucial for 
brown adipocytes to obtain glucose for immediate utilization or 
storage. We determined the messenger RNA levels of glucose 
transporters in the brown adipose tissue of rats exposed to 4°C and 
22°C for 4 hours. As expected, glucose transporter 4 (Glut4/ 
Slc2a4) was much more abundant than glucose transporter 1 
(Glutl/Slc2al) (Table 3). Glutl messenger RNA was up-regulated 
by cold exposure, p<0.001. Glut4 messenger RNA was not 
significandy changed by cold exposure as compared to the control. 

Cold-induced insulin signaling changes at the protein 
level 

We measured the levels of Insr, Irs2, Irsl, Akt-p, Akt-t, Glut4, 
and P-actin proteins in the BAT lysates from the rats exposed to 
the three conditions described in Methods (Figure 2. A and B), 
Cold and insulin-stimulated BAT lysates appeared to have more 
Insr, Irs2, Irsl, Akt-p, Akt-t, and Glut4 than controls. Irs2, Irsl, 
and Akt-p in cold group and Irsl and Akt-p in insulin-stimulated 
group were higher than controls (p<0.05). 

We explored the insulin signaling pathway in the skeletal 
muscles from the same rats which provided BAT. Since the P-actin 
is absent in skeletal muscles, we normalized the protein levels with 
total actin. Cold and insulin-stimulated skeletal muscle lysates had 
more Akt-p than controls (p<0.05) (Figure 3). The protein levels of 
Insr, Irsl, Akt-t and Glut4 were similar to the controls without 
significant changes. Irs2 signals in the skeletal muslces were very 
weak in all of conditions. 

Discussion 

In the present study, we found that 4-hours of cold exposure 
regulated multiple transcripts involved in the insulin, PKA, PI3K/ 
AKT and ERK/MAPK pathways in the BAT of rats. These four 
canonical pathways were the top pathways responding to the 
4 hours of cold exposure based on their significances compared to 
controls. It has been documented that PKA, PI3K/AKT and 
ERK/MAPK signaling pathways were altered in BAT following 
adrenergic agonist stimulation [19,20,21]. Given the fact that the 
effect of cold exposure on the BAT is substantially mediated by the 
sympathetic system [1], it is not unexpected to see the changes of 
these three canonical pathways in the BAT of rats after the 
4 hours of cold exposure. 

Our data showed that the insulin signaling pathway was the 
most affected canonical pathway responding to 4 hours of cold 
exposure. Since the up-regulated transcripts were much greater 
than the down-regulated transcripts, this change suggested that the 
direction of the insulin signaling pathway was likely up. This was 
further supported by our western blotting data of insulin signaling 
and is consistent with the chronic cold exposure study of 
Gasparetti et al [13]. 

Physiologically, it is likely that the functioning brown adipose 
tissue in normal mammals responds to cold environments by 
increasing energy consumption for the short term and increasing 
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Table 2. Predicted activation state of carbohydrate metabolism after 4 hours of cold exposure. 







Functions 


p-Value 


Predicted Activation State 


Activation z-score 


# Molecules 


metabolism of carbohydrate 


5.16E-10 


Increased 


2.643 


136 


synthesis of carbohydrate 


7.60E-10 


Increased 


2.556 


101 



doi:l 0.1 371/journal.pone.0099772.t002 



the number of mature brown adipocytes for the long term[l]. 
Carbohydrates are a major energy source for a physiological 
response to the cold in BAT [22]. 

Cameron and Smith reported that cold induced the multUocular 
cells of BAT to lose their lipid vacuoles and led to a decrease in 
their size in the first 6 to 12 hours of cold exposure, but these 
features are restored to normal by 24 hours of cold exposure. Cell 
proliferation, as estimated by the DNA synthetic index method 
(using tritiated thymidine autoradiography), appeared in the 
brown fat at 1 day of cold exposure, became maximal at 4 days 
of cold exposure, and returned to the control level by 1 6 days of 
cold exposure [23]. 

Rat studies have demonstrated that cold exposure decreased 
insulin secretion [9,10,11,13]. Improving the insulin receptor 
signaling pathway is likely an important mechanism for brown 
adipose tissue to obtain glucose more efficiently in the setting of 
cold-induced hypoinsulinemia [13]. The current transcriptome 
data matched this scenario with concurrently up-regulated crucial 
components in the insulin receptor signaling pathway including 
the insulin receptor, insulin receptor substrate I, and insulin 
receptor substrate 2. Western blotting data demonstrated 
responses which supported our data that 4 hours of cold exposure 
significantiy induced more Irs2 and Irsl than controls. These data 
were also similar to the effects of insulin on the BAT of rats kept at 
the room temperature. Even though Gasparetti et al did not 
record a significant increase in the protein level of the insulin 
receptor and Irsl, they found that the protein level of Irs2 
increased significantly. 

Our data (4°C for 4 hours) suggested that the balance between 
kinase (Insr upregulation) and protein tyrosine phosphatases (Lar 
and Ftp lb) transcripts favored phosphorylation. This correspond- 
ed to Gasparetti et al findings that 8 days of cold exposure 
significantiy enhanced the effects of insulin on phosphorylation of 
the insulin receptor, insulin receptor substrate 1, and insulin 
receptor substrate 2 in the brown adipose tissue of male rats [13]. 

The PI3K branch is well- known for regulating glucose 
metabolism in adipose tissue while insulin binds to its receptors 
[24]. Akt is activated by phosphorylation via Irs2/PI3K and Irsl/ 
PI3k in adipose tissue [24]. Our data showed that the cold 
condition significantly upregulated Pdkl mRNA level, and 
enhanced Akt-p at the protein level. As compared to the positive 
control of activating insulin signaling by administrating a dose of 
insulin to rats, we found that 4 hours of cold stimulation had 



similar effects as insulin on the BAT of rats. This finding supported 
cold stimulating the insulin signaling pathway similar to the effect 
of this insuhn dose on BAT. 

We demonstrated that the current cold condition up-regulated 
the MAPK branch. The underlying biology is most likely related 
to a marked reduction in the lipid storage in the brown adipose 
tissue as shown in our previous study after an acute cold expose 
[14]. Other data support this [23,25,26]. After lipid depletion the 
brown adipose tissue has to depend on exogenous energy sources 
such as glucose and fatty acids to maintain the thermogenesis [22]. 

Our data also demonstrated that cold exposure up-regulated the 
multiple transcripts of carbohydrate metabolism. These changes 
are consistent with prior studies showing the cold induced an 
increase in glucose uptake in brown adipose tissue [3,4,6,9,1 1,12]. 

Shivering is a common phenomenon in mice, particularly when 
first moving them from room temperature to an extreme cold 
environment such as 4°C [27]. Mizuma et al reported that a cold 
ambient temperature significantly increased FDG uptake in 
muscles and BAT of mice (18°C vs. 38°C) due to shivering [28]. 
Using a FDG micro-PET scanner, Wang et al demonstrated that 
intense FDG activity accumulated in their BAT after the 
experimental mice had been put in a cold room (4°C) for 5 hours 
(FDG was injected after 4 hours and the mice sacrificed one hour 
later) but no obvious FDG activity was identified in muscles [29]. 

In our previous experiments, 4 hours of cold exposure (4°C) did 
not significantly enhance the FDG uptake in the muscles of rats 
one-hour post FDG injection as compared to controls (p>0.05) 
[14]. This suggested that the 4°C experimental condition might 
have been easily adapted by rats or mice through a short-period of 
shivering to non-shivering thermogenesis or no shivering thermo- 
genesis at all. The Akt-p levels in the skeletal muscles of the rats 
stimulated by cold and insulin were elevated, and suggested that 
insulin signaling likely was occurring in the skeletal muscles. The 
vmderling mechanism of the cold exposure to increase the insulin 
signaling pathway in BAT and muscles should be further explored. 

Glucose transporters, often the rate-limiting step for glucose 
clearance, presented a difierent response to cold exposure. In our 
short term cold exposure, glucose transporter 4 messenger RNA or 
protein, the predominant subtype in brown adipose tissue, was not 
increased as was reported in the long term cold exposure 
experiments [13,30,31]. No significant Glut4 increase in the 
BAT and muscle lysates from cold and insulin stimulation was 
detected by western blotting. Based on increased Akt-p by cold 



Table 3. Cold-induced changes of the glucose transporters. 





Symbol 


Gene Title 


Cold (log2) 


Control (log2) 


Fold-change (log2) 


p-value 


Affymetrix 


Slc2al 


solute carrier family 2 (facilitated glucose 
transporter), member 1 


7.69 


6.85 


0.84 


0.001 


1 370848_at 


Slc2a4 


solute carrier family 2 (facilitated glucose 
transporter), member 4 


12.26 


12.19 


0.07 


0.626 


1 367989_at 



doi:l 0.1 371 /journal.pone.0099772.t003 



PLOS ONE I www.plosone.org 



June 2014 | Volume 9 | Issue 6 | e99772 



Cold Induced Response of Insulin Signaling of BAT 




RT 



Cold 



RT+ins 



(B) 



g> S 150 



9, o 



1 



i 



Glut4 
Act in 



□ Control 

□ Cold 

□ RT+lns 



#1 ntli 



InsulinR 



Akt-p 



Akt-t 



Glut4 



mechanism of regulating glucose transporter 1 in brown adipose 
tissue may be related to the overexpression of the Ras family 
[32,33]or via pS adrenergic receptor pathways [34]. 

Using primary brown adipocytes, DaUner et al reported that 
mature brown adipocytes had more Glut4 and less Glutl than 
premature brown adipocytes [34]. Activating P3 adrenergic 
receptors with norepinephrine increased Glutl (peak at 2 hours 
post activation) and decreased Glut4 expression (nadir at 5 hours 
post activation) [34]. This is similar to our cold-induced Glutl 
change in the BAT of rats, suggesting possible P3 adrenergic 
receptor involvement of Glutl upregulation. DaUner et al also 
found that Actinomycin D abolished the late phase of glucose 
uptake, but it did not abolish the early phase of glucose uptake in 
brown adipocytes stimulated by P3 adrenergic receptors and it did 
not affect insulin-induced glucose uptake in brown adipocytes 
[34]. 

Gasparetti et al found that the effects of cold on insulin signaling 
and glucose uptake in the BAT of rats were mediated by P3- 
dependent and -independent mechanisms [35]. They demon- 
strated that cold stimulation, BRL37344 (a P3 agonist) and 
SR59230A (a P3 antagonist) increased Glut4 protein levels in the 
BAT of rats [35] . However, there was a counteractive effect of the 
combination of cold stimulation and BRL37344 on Glut4 
expression, which suggested that the cold condition desensitized 
the BRL37344 signal transduction in the BAT of rats [35]. 
Gasparetti's group did not assess Glutl under their conditions. 
Further evaluation of glucose transporter efficacy under different 
conditions might help to illustrate the mechanism of cold-induced 
glucose utilization in the BAT. 

Conclusion 



Figure 3. Expression of tlie insulin receptor, Irs1, Akt-p, Alct-t, 
Glut4, and Actin in the siceietai muscles of rats with/without 
stimulation. Nine rats were fasted overnight. The next morning, they 
were randomly exposed to one of three conditions: room temperature 
+ saline (RT, control), cold temperature + saline (4X for 4 hours, Cold), 
and room temperature with insulin administration (0.75 units/kg, 
RT+lns). The preparation of BAT lysates was performed as described 
in Materials and Methods. Panel A: Protein levels of Insulin receptor, 
Irsi, Akt-p, Akt-t, Glut4, and Actin. Panel B: Columns represent mean ± 
SEM protein levels determined from "A" as a percentage control values 
of three experiments (t-test with a = 0.05, * p<0.05). The insulin 
signaling pathway in the muscles of rats post the acute cold exposure 
was likely more active than controls, as supported by increased Akt-p. 
doi:1 0.1 371 /journal.pone.0099772.g003 

and insulin stimulation, Glut4 translocation had a high likelihood 
of occurrence. We observed significantiy up-regulated glucose 
transporter 1 (Glutl), consistent with a previous report [14]. The 



This study demonstrated that 4-hours of cold exposure 
stimulated the insulin signaling pathways in the BAT of rats. 
Further exploration of the insulin signaling in the BAT of diabetic 
and obese models under cold regimens wiU provide more 
information for potential clinical translation. 
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